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Spintronics has attracted wide attention by promising novel functionalities derived from both 
the electron charge and spin. While branching into new areas and creating new themes over the 
past years, the principal goals remain the spin and magnetic control of the electrical properties — 
essentially the I-V characteristics — and vice versa. There are great challenges ahead to meet these 
goals. One challenge is to find niche applications for ferromagnetic semiconductors, such as GaMnAs. 
Another is to develop further the science of hybrid ferromagnetic metal/semiconductor heterostruc- 
tures, as alternatives to all-semiconductor room temperature spintronics. Here we present our repre- 
sentative recent efforts to address such challenges. We show how to make a digital magnetoresistor 
by combining two magnetic resonant diodes, or how introducing ferromagnetic semiconductors as 
active regions in resonant tunneling diodes leads to novel effects of digital magnetoresistance and 
of magnetoelectric current oscillations. We also discuss the phenomenon of tunneling anisotropic 
magnetoresistance in Fe/GaAs junctions by introducing the concept of the spin-orbit coupling field, 
as an analog of such fields in all-semiconductor junctions. Finally, we look at fundamental electronic 
and optical properties of GaMnAs by employing reasonable tight-binding models to study disorder 
effects. 

INTRODUCTION 

Spintronics is an interdisciplinary field driven by goals of technological innovations and fundamental discoveries. 
The two aspects, applied and fundamental, are intimately intertwined, pushing spintronics research forward. Thus 
far, the only important technological achievements of spintronics remain the giant magnetoresistance and the tunen- 
ling magnetoresistance, both based on metallic heterostructures Jj. While semiconductor spintronics has attracted 
great deal of attention, most of its device structures are still in the proposal stage, with only few device principles 
experimentally demonstrated Moving ahead with realizing potentially useful spintronic devices is a major task. 

Spintronics has also generated a great deal of fundamental knowledge. We now have firm understanding of the 
principles of spin injection, spin transport, and spin relaxation; we have learned how to manipulate spins from bulk to 
nanoscales, or how to bring ferromagnetism into useful semiconductors [3]. Important, we have learned the limitations 
of the field, sometimes fundamental in nature, sometimes technology limited. 

This article illustrates spintronics perspectives by discussing selected topics, from quantum spintronic devices, 
through spin-orbit coupling in heterostructures and its implications for tunneling and magnetocrystalline anisotropy, 
to physical properties of ferromagnetic GaMnAs. There are two important questions that these topics address in their 
specific forms: (i) What potential spintronic applications can we think of for ferromagnetic semiconductors in general, 
and for GaMnAs in particular? and (ii) Are hybrid ferromagnetic metal/semiconductor systems a viable alternative 
to all-semiconductor spintronics? The need to answer the second question stems from the limitations of the physical 
properties of useful ferromagnetic semiconductors, such as GaMnAs, whose ferromagnetism at room temperature 
remains elusive. 

Where do ferromagnetic semiconductors fit? In most proposals for spintronic devices ferromagnetic semiconductors 
provide the exchange splitting to differentiate (mainly transport) properties of spin up and spin down electrons. Here 
we discuss our proposal to take the next step: couple ferromagnetism dynamically with transport and charging, in 
resonant diodes. We discuss how such a coupling, in the quantum well made of a ferromagnetic semiconductor, leads 
to robust microwave current oscillations, at nominally dc biases. 

One great outcome of spintronics is a good understanding of spin-orbit coupling in semiconductor interfaces. Due 
to the current limitations on the Curie temperature in GaMnAs, we need to consider conventional sources of spin- 
polarization: ferromagnetic metals. In particular, ferromagnetic metal/semiconductor interfaces are attractive systems 
for controlling orbital effects by spins. We discuss here one particular interface, that of Fe/GaAs, and show that 
tunneling through such junctions leads to the phenomenon of the tunneling magnetoresistance (the tunneling resistance 
depends on the orientation of the magnetization in the plane perpendicular to the tunneling), which can be explained 
by the concept of the spin-orbit coupling field at such metal/semiconductor interfaces. Both model and first-principles 
calculations allow to explore the anisotropics stemming from this interface induced spin-orbit coupling. We also point 
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out that while the transport and magnetocrystalline symmetries both reflect the interface structure, the observed 
corresponding anisotropics have different electronic origin. 

While finding niche applications for GaMnAs is a nontrivial task, the material itself is rather fascinating. It is a 
highly p-doped semiconductor whose ferromagnetic order is mediated by holes. There are many fundamental questions 
still open: the Fermi level lies in the valence or in the impurity band? What is the effective mass of the carriers at 
the Fermi level and what is the mode of the transport? Are the states at the Fermi level fully extended or partially 
localizes? What is the value of the band gap? and so on. We attempt to answer some of those questions here with two 
reasonable tight-binding models, finding that each model explains certain experimental features in a complementary 
way, but neither provides a complete picture; a realistic model of GaMnAs is yet to be presented. 

We first introduce magnetic resonant diodes, then discuss tunneling anisotropic magnetoresistance and the spin-orbit 
coupling effects at ferromagnetic metal/semiconductor interfaces, and finally present our results and their interpreta- 
tions of our tight-binding calculations of large GaMnAs model systems. 



MAGNETIC RESONANT TUNNELING DEVICES 



Spintronic devices aim at enhancing the functionality of the existing electronic technology by exploiting the spin- 
dependent properties of solid-state systems. An efficient spin control of the charge current would allow for the 
realization of ultimate magnetoelectric devices, in which information can be processed based on the spin-state of the 
system [ij. In this context, band-engineered magnetic semiconductor heterostructures offer promising perspectives, 
since they provide opportunities for controlling and tuning spin-dependent magneto-transport properties In 
particular the rapid development of novel dilute magnetic semiconductors (DMS) in the last decades [1, HTSl' 
which semiconductors are made magnetic by doping with transition metal elements, has considerably enlarged the 
class of suitable materials for the growing of heterojunction systems. 

Magnetic heterostructures typically consist of stacked layers of only a few nanometers of both magnetic and non- 
magnetic semiconductors. Along the growth direction the system size becomes comparable to the phase coherence 
length of the propagating carriers, allowing quantum interference effects to govern the vertical transport properties. 
According to the band-mismatch of the different semiconductor materials, some of the layers constitute energetic bar- 
riers for the incident carriers, which could only be penetrated by tunneling. Nevertheless, in double and multi-barrier 
systems transmission probabilities up to nearly 100% can be found for certain incident carriers energies, which actu- 
ally correspond to the discrete electronic spectra of the formed quantum wells; an astonishing effect, termed resonant 
tunneling. 

From a generic point of view, resonant tunneling structures act like an energy filter, allowing only carriers with 
resonant energies to pass through. The transport becomes strongly sensitive to the relative alignment of the electronic 
spectra in the leads and the wells, which can be modified by external bias voltages. Even small energy splittings for 
different spin states in the quantum wells, e.g., induced by an external magnetic field or an exchange field, can result in 
strong ramifications in the transport characteristics, as has been observed experimentally in paramagnetic ^ as well 
as ferromagnetic [7,, ^8^] DMSs quantum well structures. Due to the spin-dependent transmission magnetic resonant 
tunneling devices have been proposed for realizing efficient spin valve, spin filtering, spin switching, and spin detecting 
devices , which form important building blocks for a spin-based information technology. 

In the following we show representative examples of capabilities of magnetic heterostructures by describing (i) the 
digital magneto resistance (DMR), which appears in serially connected resonant tunneling diodes (RTDs), and (ii) the 
dynamic interplay of transport and magnetic effects in ferromagnetic RTDs, which can lead to selfsustained current 
oscillations at a nominally steady bias voltage. 

A nice example, which illustrates that the current characteristics can be strongly modified by changing the magnetic 
properties of resonant tunneling structures, is given by a double barrier structures comprising a ferromagnetic quantum 
well and a ferromagnetic emitter, as illustrated in Fig. [Ijb). Assuming that the magnetization of the quantum well is 
"soft" the relative alignment of the emitter and well magnetization directions can be changed by external magnetic 
fields. Selfconsistent numerical simulations for GaMnN-like systems ll| reveal a strong dependence of the current- 
voltage (IV) characteristics on the relative magnetization angle, as shown in the lower panel of Fig.[TJ In particular, 
the peak current is considerably decreased if the the well magnetization is tilted away from the parallel alignment 



with the emitter's magnetization. This effect can be used to realize the proposed DMR-effect [11|, |12| in so-called 
monostable-bistable logic elements (MOBILEs) [H, [11] . 

A conventional MOBILE consists of two serially connected RTDs, a load and a driver, as shown in Fig. [Ija). The 
principle of the device operation relies on the nonlinear IV-characteristics of RTDs. By drawing the circuit load line 
diagram a single monostable output voltage for low input voltages is obtained. However, at high input voltages a 
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FIG. 1: Top panel: (a) circuit diagram of a magnetic MOBILE. The load is a conventional RTD, whereas the driver is a RTD 
with a ferromagnetic emitter and quantum well (b). The well magnetization can be tilted by an external magnetic field B by the 
angle (p. Lower panel: Selfconsistent current-voltage characteristic of the ferromagnetic driver-RTD for different orientations of 
the quantum well magnetization, indicated by the angle 0. 



bistable working point regime appears. Which of the two working points is realized actually depends on the difference 
of the load and driver peak currents: if the load peak current is smaller (higher) than the driver's peak current than 
the output voltage is low (high). In a nutshell: a MOBILE can compare the peak currents of the load and driver 
RTDs, yielding a digital binary output signal. 

In a magnetic variant of the MOBILE the driver is replaced by a magnetic RTD, e.g, of the kind as discussed 
above. If the quantum well magnetization is tilted by an external magnetic field the magnetic driver's peak current is 
continuously decreased. At some critical angle the driver's peak current becomes smaller than the load peak current, 
which results in a discrete jump of the output voltage from low to high. This effectively realizes a digital magneto 
resistance: the total device resistance makes a discontinuous step by continuously changing the magnetic properties 
of the structure, i.e, in this case the tilting angle of the quantum well magnetization. 

The DMR-like effect has, in fact, already been experimentally demonstrated in a somewhat different device config- 
uration, in which a metallic giant magneto resistance (GMR) element is shunted to a nonmagnetic driver diode [l5|. 
In terms of applications DMR-devices could be used as very fast read heads of hard disks, in which the magnetically 
stored information is directly converted into a binary electrical signal. Furthermore, ferromagnetic RTDs can have 
the advantage to be nonvolatile upon the loss of power, since the state of the device is stored in the magnetization 
direction of particular layers of the structure. Such nonvolatile devices are attractive for fast random access memory 
applications or for reprogrammable logics, in which the logical circuit function is altered by changing the magnetic 
state of the device. 

In addition to magnetically tunable steady state properties, interesting dynamic phenomena can occur in magnetic 
heterostructures according to a peculiar feedback process, in which the magnetic ordering of the quantum well affects 
back the resonant tunneling current. In DMSs the ferromagnetic order is mediated by the itinerant carriers, which 
gives rise to a close interdependence of the transport and magnetic properties. Indeed, several experiments succeeded 
in generating ferromagnetism by tailoring the carrier density by electrical or optical means 1^, 17 1. 

In resonant tunneling structures the particle density in the well strongly depends on the resonant tunneling condi- 
tions, which implies that the critical temperature of magnetic quantum wells should be tuneable by external biases, as 
theoretically predicted in the framework of a mean field model of the well magnetism 18, 1^. These studies revealed 
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FIG. 2: Top panel: transients of the oscillating spin-resolved current density at a fixed dc-bias. Lower panel: qualitative 
explanation for the occurrence of self-sustained oscillations, (a) The in-tunneling spin-up carriers increase the spin polarization 
5, which causes an increasing well splitting A. (b) When the spin up level crosses the cutoff energy the total particle number 
n and consequently the electrostatic potential Uw decreases, (c) This pushes the spin-up level even deeper into the cut-off 
region leading to a fast decrease of the spin polarization and, hence, of the well splitting, bringing the spin up level back to the 
emitter's supply region (d), where the whole process restarts again. 



that the exchange splitting of the carriers subbands essentially depends (i) on the spin polarization of the particle 
density in the well and (ii) on the overlap of the subband wave function with the magnetic impurity density profile. 

Based on these findings we recently proposed that the combined nonlinear feedback of Coulomb and magnetic 
interaction can lead to self-sustained high-frequency oscillations of the tunneling current and the quantum well mag- 
netization over a large window of nominally dc-bias voltages [2Q[, as illustrated in the upper panel of Fig. [2l The 
occurrence of these magnetoelectric oscillations needs a built-in energy cut-off Scut of the emitter tunneling rate. 
This might be realized, e.g, by a cascaded left barrier, which exponentially suppresses the tunneling of carriers with 
energies smaller than i?cut- 

A qualitative discussion and explanation for the appearance of the self-sustained oscillations is given in the lower 
panel of Fig. [2| the electrostatic feedback behaves like an "inertia", allowing the spin level with energy Eg to get 
deeper into the region, in which no carriers can tunnel from the emitter side anymore {Eg < E^ut), whereas the 
magnetic feedback is required to push the spin level back into the emitter's supply region {Eg > E^ut)- The frequency 
of the oscillations can be modified by the applied dc-bias voltage, which suggest applications of ferromagnetic RTDs 
as tunable high-power oscillators. 

These two examples of magnetic resonant tunneling structures, mainly chosen on the basis of the authors personal 
research interest, illustrate that magnetic heterostructures have a promising potential for realizing fully integrated 
magnetoelectronic devices as well as for gaining a deeper understanding of the mechanism of carrier-mediated ferro- 
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FIG. 3: Schematics of a TAMR MTJ with a single magnetic layer (the top layer). The vector n points in the direction of the 
magnetization and [x] denotes the reference axis for 0. 

magnetism. Especially the possibilities of ferromagnetic multi-barrier systems, which naturally provide a rich variety 
of compositional configurations, are still largely unexplored. 

TUNNELING ANISOTROPIC MAGNETORESISTANCE IN MAGNETIC JUNCTIONS 

Spin valve devices are layered heterostructures combining ferro- and paramagnetic layers. Such devices are at the 
heart of the spintronics and are important for the investigation and measurement of spin injection, spin relaxation, 
and spin polarization. The spin valve devices find also multiple applications in the design of read head sensors and 
data storage devices, such as magnetic random access memory. The operation of conventional spin valve devices relies 
on the dependence of the magnetoresistance on the relative orientation of the magnetization in the ferromagnetic 
leads. The manifestation of such a phenomenon in tri-layer systems composed of two magnetic layers separated 
by a non-magnetic metallic spacer is called the giant magnetoresistance effect, while in magnetic tunnel junctions 
(MTJs) composed of an insulating or semiconducting barrier sandwiched between two magnetic electrodes is referred 
to as the tunneling magnetoresistance (TMR) effect. It has been observed however, that the magnetoresistance 
in MTJs may also depend on the absolute orientation of the magnetization in the ferroma gnet ic leads 2l|, 



. The tunneling 
. These authors 



24l |. This phenomenon is called the tunneling anisotropic magnetoresistance (TAMR) effect [21 
magnetoresistance in GaMnAs/GaAlAs/GaMnAs tunnel junctions was theoretically investigated 
predicted that, as a result of the strong spin-orbit interaction, the tunneling magnetoresistance depends on the angle 
between the current flow direction and the orientation of the electrode magnetization. Thus, a difference between the 
tunneling magnetoresistances in the in-plane (i.e., magnetization in the plane of the magnetic layers) and out of plane 
configurations of up to 6% was predicted for large values of the electrode spin polarization [23| . Here we refer to this 
phenomenon as the out-of-plane TAMR (see Fig. [3]). Recent ab initio calculations in Fe/MgO/Fe magnetic tunnel 
junctions (MTJs) predict an out-of-plane TAMR ratio of about 44 % [25]. On the other hand, the in-plane TAMR 
effect refers to the changes in the magnetoresistance when the in-plane magnetization of the ferromagnetic layer(s) is 
rotated around the direction of the current flow (see Fig. [3]). 

It is remarkable that the TAMR is present even in MTJs in which only one of the electrodes is magnetic and the 
conventional TMR is absent. Thus, a new generation of TAMR-based devices could have some advantages with respect 
to the devices whose operation is based on the TMR effect. In contrast to the conventional TMR-based devices, which 
require two magnetic layers for their operation, TAMR-based devices can, in principle, operate with a single magnetic 
lead. Thus the TAMR effect could open new possibilities and functionalities for the operation of spintronic devices. 
The TAMR effect may have an influence on the spin-injection from a ferromagnet into a non-magnetic semiconductor. 
Therefore, in order to correctly interpret the results of spin-injection experiments in a spin-valve configuration, it is 
essential to understand the nature, properties, and origin of the TAMR effect. 

The first experimental observation of the TAMR effect was performed in (Ga,Mn)As/A10x/Au heterojunctions, 
where an in-plane TAMR ratio of about 2.7% was found [2lj. Experimental investigations of the in-plane TAMR in 
(Ga,Mn)As/GaAs/(Ga,Mn)As :22] and (Ga,Mn)As/ZnSe/(Ga,Mn)As fSs*] tunnel junctions in which both electrodes 
are ferromagnetic have also been reported. In the case of (Ga,Mn)As/ZnSe/(Ga,Mn)As the in-plane TAMR ratio 
was found to decrease with increasing temperature, from about 10% at 2 K to 8.5% at 20 K [2^. This temperature 
dependence of the in-plane TAMR is more dramatic in the case of (Ga,Mn)As/GaAs/(Ga,Mn)As, for which a TAMR 
ratio of order of a few hundred percent at 4 K was amplified to 150 000% at 1.7 K j22|]. This huge amplification of 
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FIG. 4: (a) Angular dependence of the TAMR in a Fe/GaAs/Au MTJ for different values of the bias voltage V^^^^. Solid lines 
correspond to theoretical results while symbols represent the experimental data, as deduced from Ref. [s?]. (b) Extracted bias 
dependence of ai. 



the in-planc TAMR was suggested to originate from the opening of an Efros-Shklovskii gap [26| at the Fermi energy 

Measurements of the TAMR in p+-(Ga,Mn)As/n+-GaAs Esaki 



when crossing the metal-insulator transition 
diode devices has also been reported 



22| 



27[ . In addition to the investigations involving vertical tunneling devices the 



TAMR has also been studied in break junctions [28 



271 130| and nanocontacts [31 1 



29|, nanoconstrictions 

In addition to the currently low Curie temperature ferromagnetic semiconductors, the TAMR has recently been 
experimentally investigated in Fe/GaAs/Au MTJs Co/AlO^r/Au MTJs [33*1, and in multilayer-(Co/Pt)/A10:r/Pt 
structures |34j . 

In what follows we shall focus our discussion on the case of the in-plane TAMR (for brevity we will refer to it just 
as the TAMR effect) recently measured in Fe/GaAs/Au MTJs [32|. In such an experiment the magnetization of the 
ferromagnetic Fe electrode was rotated in- plane by an angle </> with respect to the [110] crystallographic direction and 
the magnetoresistance R{4>) = R{9Q°, (f>) measured as a function of (see Fig. [3]). The strength of the TAMR is given 

by i 



TAMR[no](0) 



i?(0) 



(1) 



The angular dependence of the TAMR in a Fe/GaAs/Au tunnel heterojunction is shown in Fig. lUJa) for different 
values of the bias voltage V^^.^^. The symbols represent the experimental data (conveniently mirrored) as deduced 



from Ref. [34]. Two fundamental features in the behavior of the TAMR can be appreciated in Fig. H^a): (i) the 
experimentally measured TAMR exhibits a two-fold symmetry with respect to the magnetization direction in the 
ferromagnet, (ii) by varying the bias voltage the TAMR can be inverted (i.e., its sign can be changed). In order to 
explain these experimental findings we have proposed a model (for details see Ref. 0) in which the two- fold symmetry 
of the TAMR in ferromagnet/semiconductor/normal metal junctions emerges from the interference of Dresselhaus and 
Bychkov-Rashba SOCs. The combined action of these two spin-orbit interactions leads to net, anisotropic SOC with 
a two-fold symmetry which is compatible with the C2V symmetry of the semiconductor interfaces and manifest in 
the tunneling magnetoresistance. The Dresselhaus SOC results from the bulk inversion asymmetry of the zinc-blende 
semiconductor (GaAs in this case) and is characterized by the Dresselhaus parameter 7 which has a finite value in 
the semiconductor region (for GaAs 7 « 24 eV Q) and vanishes elsewhere. 

For the small voltages considered in the experiment the main contribution to the Bychkov-Rashba SOC comes from 
the strong potential gradients at the interfaces rather than from an external electric field. The interface Bychkov- 
Rashba SOC arises due to the semiconductor interfaces inversion asymmetry Q and is characterized by the parameters 
ai and which measure the strengths of the SOC at the left (ferromagnet/semiconductors) and right (semiconduc- 
tor/normal metal) interfaces, respectively. These parameters are not known for metal/semiconductor interfaces and 
have to be extracted from available experimental data j32l | or from first-principles calculations (see the following 
section). For the case of an Fe/GaAs/Au MTJ it turns out that the size of the TAMR is dominated by the parameter 
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ai corresponding to the ferromagnet /semiconductor interface j2l ■ Since the values of ar are irrelevant for the TAMR, 
one can set this parameter, without loss of generality, to zero. This leaves ai as a single fitting parameter when 
comparing to experiment. The theoretical results are shown by solid lines in Fig. |31 The agreement between theory 
and experiment is very satisfactory. The values of the phenomenological parameter ai are determined by fitting the 
theory to the experimental value of the TAMR at ^ = 90°. This is enough for the theoretical model to reproduce the 
complete angular dependence of the TAMR, demonstrating the robustness of the model. By performing the fitting 
procedure for different values of the bias voltage V^^.^^ the bias dependence of a/ can be extracted. The results are 
shown in Fig. [UJb). The obtained values of a; have reasonable magnitudes, not too different from known values in 
semiconductor interfaces Q. For a bias voltage of about 50 mV the parameter a; « and the TAMR vanishes 
[compare Figs.lU^a) and (b)]. At a bias slightly bellow 50 mV the interface Bychkov-Rashba parameter changes sign, 
resulting in the inversion (change of sign) of the TAMR, as seen from Fig.UJa). 

In order to understand the main qualitative behavior of the TAMR, a simplified phenomenological model was 
proposed in Refs. 0, [l^l- Such a model is based on a perturbative expansion of the transmissivity in powers of n.w, 
with n denoting the unit vector defining the magnetization direction and w = {—diky + 'ykx, dikx — jky, 0) being the 
effective spin-orbit coupling field. Here a and 7 represent the average values of the Bychkov-Rashba and linearized 
Dresselhaus parameters, respectively (note that ai oc ai and 7 ex 7). It was then shown, on the basis of very general 
symmetry considerations, that up to second order in the strength of the SOC field the angular dependence of the 
TAMR is determined by the relation (for details see Ref. [1]) 

TAMR[no] « aa[cofi{2cb) - 1], (2) 

which is consistent with the experimental data and the full theoretical calculations [see Fig. Hl^a)]. One can see from 
Eq. ^ that bias-induced changes of the sign of ai lead to the inversion of the TAMR. It is also clear that when 
Q!;7 = 0, the TAMR is suppressed. Therefore, for a finite in-plane TAMR both Bychkov-Rashba and Dresselhaus 
SOCs must necessarily be present. This is in contrast to the case of the TAMR in the out-of-plane configuration for 
which the presence of only Bychkov-Rashba coupling suffices 13, 1^ ■ 



AB-INITIO STUDY OF SPIN-ORBIT COUPLING EFFECTS IN FE/GAAS INTERFACES 

The spin-orbit interaction couples electron spin to the orbital moment. In magnetic systems the spin-orbit inter- 
action is the primary source of the magnetocrystalline anisotropy. In Fe films grown on semiconductors with a cubic 
zinc-blende structure, e.g., ZnSe, GaAs or InAs, a uniaxial magnetic anisotropy has been reported [36| . This comes 
rather as a surprise since both the Fe and the semiconductors have cubic symmetry. The easy axis has been found 
along [110] or [110] directions. The origin of the different easy axis orientations has been attributed to the different 



semiconductor surface reconstructions. Kneedler et al. [37| showed that strain in the Fe film grown on GaAs(OOl) 



surface has negligible effect on the easy axis orientation indicating that it is independent of the surface reconstruction. 
Moosbiihler et al. ^^8\ investigated (4 x 2) and (2 x 6) surface reconstructions and pointed out that easy axis is always 
along [110] with practically the same amplitude for different Fe film thicknesses. It is of importance to introduce a 
realistic model for the origin of the uniaxial magnetic anisotropy as well as to relate it to the TAMR effect, whose bias 
induced reversal is well described, on the one hand, by the interference between Bychkov-Rashba and Dresselhaus 
spin-orbit couplings 0, [s^] and, on the other hand, by the reversal of the spin polarization [391 ] . The interface structure 
also affects the Schottky barrier height [13] and consequently spin injection to the GaAs. 

As the lattice constant of the bulk GaAs (5.65 A) is almost twice the lattice constant of the bcc Fe (2.87 A), the 



smooth epitaxial growth of Fe on a GaAs(OOl) surface is possible [3a,|4l|. Erwin et al. [42| studied, using density 
functional theory, the stability of ideal 1x1 Fe/GaAs interfaces. Their study revealed that the thickness of the Fe 
film as well as the type of the GaAs surface termination (As- or Ga-like) affects the interface stability. Three types 
of 1 X 1 Fe/GaAs interface models have been investigated [i^: (i) model A - flat interface [Figl5ja)]; (ii) model B 
- partially intermixed interface [FiglUb)]; (iii) model C - fully intermixed interface [Figl5l(c)]. When more than 
two atomic layers of Fe are deposited on a GaAs(OOl) surface, the flat (model A) and partially intermixed (model B) 
interface models are more stable than the fully intermixed ones (model C). As the pd hybridization lowers the cohesive 
energy [i^], the antibonding bonds between Fe and As are more stable than the Fe-Ga bonds. Interestingly, while the 
density functional calculations imply that the As-terminated flat (model A) interface is more stable then the partially 
intermixed (model B) (4^ one, a recent Z-contrast scanning transmission electron microscopy experiment reported a 



single plane of alternating Fe and As atoms (model B) in an Fe/AlGaAs interface [4J]. 

In the following we discuss electronic properties of an ideal slab structure of the flat interface model (model A) 
between Fe and As-terminated GaAs(OOl). The slab encased in vacuum of the 6A thickness contains 9 atomic planes 
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FIG. 5; Schematics of 1 x 1 Fe/GaAs interface models (a) model A, (b) model B and (c) model C. (d) difference of d-projected 
local density of states of two nonequivalent Fe atoms at the flat interface (model A) for magnetization along [110]. 



of GaAs and three atomic planes of bee Fe on top of the GaAs(OOl) surface. The orthorombic symmetry of the system 
is reduced due to the spin-orbit coupUng, depending on the orientation of the magnetization in the Fe. The presence 
of the much debated uniaxial magnetic anisotropy can be explained in this flat interface model by means of directional 
covalent bonds at the interface [45|. The interface breaks the natural continuation of the sp^ tetrahedral bonds of 
the zinc-blende structure. Because of that, in what follows two inequivalent groups of Fe atoms at the interface can 
be distinguished. The first group contains Fe atoms Fel [see FiglSJa)], which are in "bonding" positions, e.g., the 
sites where Ga atoms would have been situated. The second group contains Fe2 atoms that correspond to voids in a 
zinc-blende structure. These nonequivalent Fe atoms show similar local densities of states (LDOS), but if one plots 
the difference of the d-projected LDOS, see FigOJd), the couphng of the Fe film to the GaAs substrate appears to 
be clear. The d-states about 0.5 eV below the Fermi level {E-p) are split in the case of Fel due to bonds between 
hybridized sp^ states (As) and d-states (Fel) to bonding (1.4 eV bellow E-p) and antibonding states (0.4 eV above 
Ep). Thus, by means of LDOS one finds the answer that uniaxial magnetocrystalline anisotropy has its origin in 
anisotropic interfacial bonds [4^, which represents a microscopic mechanism causing the magnetic anisotropy in the 
Fe/GaAs heterostructure. It has been shown in case of Zn-terminated Fe/ZnSe(001) system that all Fe atoms in 
"bonding" positions contribute to uniaxial magnetic anisotropy with easy axis along [110] direction while atoms in 
void- like positions along [110] direction with 40-times lower magnitude [45|. In the case of Se-terminated interface the 
easy axis points along [110]. In the case of Fe/GaAs heterostructure one assumes similar trends since the splitting of 
the bonding states is analogous. 

In addition to the magnetic anisotropy, a novel transport phenomenon has been discovered in Fe/GaAs heterostruc- 
tures: a tunneling anisotropic magnetoresistance (TAMR) effect, as discussed in the previous section. This effect, 
which has been observed in Fe/GaAs/Au junctions [s^] with epitaxial Fe/GaAs interfaces, has challenged our under- 
standing of ferromagnetic metal/semiconductor junctions. Indeed, a typical spin- valve like behavior has been observed 
for a junction that contains only a single ferromagnetic layer; at least two layers are needed for conventional spin 
valve phenomena. To investigate the origin of the TAMR effect have proposed a simplified model of a ferromag- 
net / semiconductor heterostructure based on general symmetry considerations 0, [s^] . The spin-orbit coupling in the 
junction is manifested through linear (in momentum) Bychkov-Rashba and Dresselhaus terms, which reflect, respec- 
tively, the structure and bulk inversion asymmetry of the Fe/GaAs heterostructure. Interference of these spin-orbit 
couplings leads in [001] quantum wells to the C2v symmetry. The same symmetry then appears in the tunneling 
magnetoresistance, which is a function of the magnetization direction 32]. In addition, we consider a spin splitting. 
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which describes the sphtting of the corresponding spin states due to the strong exchange sphtting in the Fe film. The 
eigenspectrum, in the hmit of a strong exchange sphtting (weU justified since the sphtting in Fe is about 1.5 eV, much 
more than the sphtting due to the spin-orbit interaction) close to the Brillouin zone center is 

<{k) = el") ^a{k K sin(0 - 0) + 7„ sm{^ + 0)]) . (3) 

Since we have allowed for band dependent Bychkov-Rashba (a„) and Dresselhaus (7„) parameters, we can identify the 
effect of the spin-orbit coupling on the parabolic band el") as a shift that depends on the spin (cr = ±1), the sign of an 
and 7„, the orientation of the magnetization {(jj) as well as the orientation of the momentum (9). Performing a fitting 
of ab-initio data to the eigenspectrum Eq.Q we have obtained effective Bychkov-Rashba and Dresselhaus parameters. 
We have found that the Bychkov-Rashba parameters change sign within the energy window of about 250 meV above 
and below E-p. Since the sign of TAMR depends on the product [see Eq.(l2])] we can conclude that the bias induced 
TAMR reversal, observed in Ref. [s^, can be explained by the change of the sign of the Bychkov-Rashba parameters. 

It would be reasonable to expect that the magnetocrystalline anisotropy can be influenced by an applied bias in 
the same way as the tunneling anisotropic magnetoresistance. As we have seen, though, the two effects originate 
from different electronic states: magnetocrystalline anisotropy is given by the states deep below the Fermi level, while 
TAMR is determined by the propagating states close to the Fermi level. The bias dependence of the TAMR means that 
states within a transport window (given by the bias) contribute. As these states have in principle different spin-orbit 
field parameters, upon averaging over the transport window we obtain different resistance at different biases; we are 
probing different states. Such an averaging would not change the magnetocrystalline anisotropy which is dominated by 
the spin-orbit effects of the deep lying states. However, one could expect effects due to the bias dependent electronic 
deformation of these states. 



TIGHT-BINDING MODELING OF ELECTRONIC AND OPTICAL PROPERTIES OF GAMNAS 



Diluted magnetic semiconductors have been subject to very active research during the past years due to their 
promising role for new spin-electronic devices Q . GaMnAs is a much studied material system among these magnetic 
semiconductors as its Curie temperature is rather high, reaching about 170 K [3, 5]. The host material, GaAs, is a 
III-V fee semiconductor with a band gap of 1.5 eV; the Mn impurity has a binding energy of 0.11 eV when it replaces 
a Ga atom, see Fig. [S) Since Mn is short one p electron compared with Ga at the Asp valence levels, substitutional 
Mn impurities act as acceptors that carry a local magnetic moment leading to a hole mediated ferromagnetism [i^ . 

Despite efforts to understand the bulk properties of ferromagnetic GaMnAs both from an experimental and theo- 
retical point of view a number of questions remain open. It is still an unsettled issue whether the free carriers (holes) 
reside in an impurity band which is detached from the host valence band or in the valence band itself. To this end 



various absorption experiments [471 148| . |49| were performed. However, the Mn concentration threshold for which the 
impurity band, that appears for very low Mn concentrations, eventually merges with the host valence band could 
not uniquely be determined. Measurements of the Fermi energy in GaMnAs and the size of the gap as a function of 





d(E) [a.u.] L 



FIG. 6: Crystal structure of GaAs with substitutional Mn impurity (left panel) and a calculated tight-binding band structure 
and the density of states (right panel). 
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FIG. 7: Density of states for Masek's model (positive sign) and Tang, Flatte's model (negative sign). The solid lines show 
the corresponding results for clean GaAs. INSET: Gap size A (triangles: Masek - upwards; Tang, Flatte - downwards) and 
Fermi energy Ef (Masek - crosses; Tang, Flatte - stars). The dashed line represents the Fermi energy of p-GaAs with the 
corresponding number of holes added. 



the impurity concentration are presented in 50, 51 1- These characteristics are important for the design and study 
of spintronic devices as well as for spin injection experiments. Further open questions concern the scattering; and 
localization properties of the holes caused by the disorder in the material. Besides the impurity band picture [49i] , an 
interpretation based on a merged impurity and valence band together with the existence of localized states in the band 
tail could also explain the absorption experiments [s^] . A wide range of theoretical models describing the electronic 
properties of GaMnAs has been reviewed in 5 L Among these are first-principles calculations 5^ 54 1 , effective single 
particle tight-binding approaches 5^ 5^ 57^5^, tight-binding approaches in combination with percolation theory 
[Hil [Hoi . [6l| . dynamical mean field theories |62{. effective theories based on fc • p Hamiltonians 63, 64, 65, and 
large-scale Monte-Carlo studies of a real space Hamiltonian [67l |. 

In our studies we address these open questions using a phenomenological multi-band tight-binding simulation 
method for ferromagnetic bulk GaMnAs. This material specific method treats the Mn impurities in a non-perturbative 
way and can be applied to system sizes of up to 2000 atoms (nowadays) with reasonable computational effort. It is 
based on well-known sp^ tight-binding models for the host GaAs material which include spin orbit coupling [gl, |6^ . 
The disorder effects are modeled by changing the on-site and hopping terms of those Ga sites that are replaced by 
the Mn impurities. We compare two recent proposals, one due to Masek and the other due to Tang and Flatte. The 
first approach is based on a second-nearest neighbor approximation for the host GaAs [6^ and determines the Mn 
tight-binding parameters by means of first principles calculations in a Mn concentration regime around 10% (ssl. [70|. 
The second approach remains within the nearest neighbor approximation [68| . It is based on a fit of two physically 
relevant tight-binding parameters to reproduce the binding energy of a single Mn impurity in the GaAs host material 
" ssl 71 1 . In either case an effective single-particle model is obtained which treats the carrier-carrier interactions in a 



mean field approximation. 

As it is not clear yet which of the two different parameter sets (if any) is better suited for the description of available 
experimental data we have investigated the density of states, the position of the Fermi energy, the size of the gap, the 
localization properties of the wavefunctions and the optical conductivity of GaMnAs for both models and compared 
the results to experimental findings [ll] . 

First, we present the results of both models for the density of states in the dependence on the Mn concentration 
X. In Fig. [7| the density d{E,x) is shown for the spin up states only as the spin down contributions follow the GaAs 
density very closely. Because the two models are based on different parameter sets for clean GaAs the resulting 
densities also show a slightly different behavior, in particular for the conduction band. 

Photoconductivity experiments [lH on GaMnAs heterostructures yield a gap A « 1.41 eV for x « 3 % with 
being close to the impurity level. Scanning tunneling microscopy measurements [s^l resulted in A « 1.23 eV for a 
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similar concentration. However, a quantitative comparison with the experimental data and our results is not possible 
due to the uncertainty in the number of Mn interstitials in the experimental data. 

For Masek's model the density of states shows only minor modifications due to the increasing disorder, see Fig. [T] 
There is no detached impurity band forming in the relevant range of Mn concentrations x. The Fermi energy decreases 
into the host valence band as the Mn impurities increase the number of holes. If the Fermi energy is measured from 
the top of the GaMnAs valence band as a function of the Mn impurity one finds a close similarity to clean p-GaAs 
with the corresponding number of holes added, see inset of Fig. [71 For the second model due to Tang and Flatte the 
situation is completely different. An impurity band, which merges with the host valence band at Mn concentrations 
> 1% can clearly be identified, see lower curves of Fig. [71 A large number of valence states are shifted into the host 
gap resulting in a Fermi energy which increases compared to the top of the host valence band. The gap eventually 
disappears at a; « 5 %; this is not consistent with experiment. 

Besides the density of states, the localization properties of the states around the Fermi energy play an important 
role for transport phenomena. We studied these properties by means of the participation ratio and found that the 
states around the Fermi energy are extended in both models. Similar analysis for energies closer to the band edge 
leads to a different result for the model by Tang and Flatte. At these higher energies in the GaMnAs valence band 
tail one finds strongly localized states . 

The optical conductivity (t(w) is closely related to the absorption and can be directly accessed in experiments [48, 
iil . TSj . However, there are still some discrepancies in the results and the interpretation of the data concerning the low 



frequency limit, the position of the intervalence band peaks, and the deduced effective masses and mean free paths. 
The two different models we studied result in two qualitatively different line shapes for Masek's model shows 

a Drude peak at zero frequency while no such peak is observed in Tang and Flatte's model. However, despite all the 
qualitative differences of the two models, i.e. concerning the existence of an impurity band, the calculated effective 
masses are similar being within the range of 0.4 - 1.0 of the free electron mass. 

We conclude that phenomenological tight-binding methods are well suited for the simulation of material properties 
of ferromagnetic semiconductors as they include the disorder in a non-perturbative way. However, neither of the two 
proposed parameterizations studied by us describes satisfactorily available experimental data. They both describe 
some properties correctly while failing at others. 
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